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ABSTRACT 



Context. In order to get a significant Zeeman signature in the polarised spectra of a magnetic star, we usually 'add' the 
contributions of numerous spectral lines; the ultimate goal is to recover the spectropolarimetric prints of the magnetic 
field in these line additions. 

Aims. Here we want to clarify the meaning of these techniques of line addition; in particular, we try to interpret the 
meaning of the 'pseudo-line' formed during this process and to find out why and how its Zeeman signature is still 
meaningful. 

Methods. We create a synthetic case of line addition and apply well tested standard solar methods routinely used in 
the research on magnetism in our nearest star. 

Results. The results are convincing and the Zeeman signatures well detected; Solar methods are found to be quite efficient 
also for stellar observations. We statistically compare line addition with least-squares deconvolution and demonstrate 
that they both give very similar results as a consequence of the special statistical properties of the weights. 
Conclusions. The Zeeman signatures are unequivocally detected in this multiline approach. We may anticipate the 
outcome that magnetic field detection is reliable well beyond the weak-field approximation. Linear polarisation in the 
spectra of solar type stars can be detected when the spectral resolution is sufflciently high. 



1. On the Multiline Zeeman Signature and 
Zeeman Doppler Imaging concepts 

In the 1980s it became clear that cool stars, mainly rapid 
rotators, exhibit solar type activity. Rapid rotation is one 
of the essential ingredients to the stellar dynamo. While 
there is not yet any completely satisfactory theory or model 
of the solar or the stellar dynamo, it is quite clear what 
drives such a process, viz. convection, rapid rotation and 
differential rotation. Tentatives to detect the Zeeman effect 
in these active stars however failed for two main reasons: 

— The magnetic configuration of the stellar magnetic fields 
can be quite complex. The simultaneous appearance of 
opposite magnetic polarities may lead to cancellation of 
the respective contributions to the integrated polarisa- 
tion signal in the stellar spectrum. 

— The broadening of the spectral lines due to even moder- 
ate stellar rotation rates can become a serious handicap 
when measuring the magnetic field of the star. 

Consideration of these two aspects was at the object of 
Semel (1989) where the term Zeeman Doppler Imaging 
(ZDI) stood for the detection of a Zeeman signature thanks 
to the Doppler effect. The latter can help to disentangle the 
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contributions from opposite magnetic polarities and their 
respective opposite polarisation which otherwise may can- 
cel. Thus, the combination of Zeeman and Doppler effects is 
one of the reasons why the detection of the Zeeman effect 
in fast rotating active solar type stars becomes possible. 
Nowadays, however, the term ZDI is often understood as 
magnetic mapping of the stellar surface, based on the inver- 
sion of Stokes profiles (frequently only / and V) sampled 
over a full period of rotation. 

In order to overcome possible misunderstandings, we 
therefore propose the term Multiline Zeeman Signature 
(henceforth MZS) to denote a method for just the detec- 
tion of a mean Zeeman signature using numerous spectral 
lines. We shall not attribute any direct physical meaning 
to the MZS. We simply require the application of an oper- 
ator (or a detector 0) to a polarised spectrum, creating 
a particular MZSq. 

Nowadays, there are several techniques that take into 
account many spectral lines to obtain some particular MZSs 
that trace stellar magnetic fields. One of the most ele- 
mentary methods is the so called line addition technique 
which consists of simply adding up many observed spec- 
tral lines (Semel 1989, Semel & Li 1996). Later on, Donati 
et al. (1997) introduced the Least Squares Deconvolution 
(LSD) method. More recently, methods based on Principal 
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Component Analysis (PCA) have been introduced by Semel 
et al. (2006) and Martinez Gonzalez et. al (2008). In this 
paper we concentrate on the proper definition of the line 
addition technique and its application to any of the Stokes 
parameters. We present the advantages of using this tech- 
nique as compared to LSD and we describe the correct way 
to infer physical properties of its MZS. 

2. The MZS through line addition technique as 
compared to the LSD profile 

2.1. Line Addition 

The first successful MZS tentatives consisted in adding up 
the circular polarisation of selected spectral lines. In the 
first efforts, Semel (1989) recommended the de-blending of 
spectral lines prior to addition, following a relatively tough 
procedure. However, Semel & Li (1996) showed that simple 
coherent addition of spectral lines gave satisfactory results. 
This is due to the fact that the the position of the blends 
along the spectral lines is not systematic and they add inco- 
herently. Each spectral line may appear several times, once 
in the coherent addition, and then each time when the line 
is "blending" in the field of another nearby line. However, 
in the latter case it appears in a non-coherent way and al- 
ways in another position (in wavelength). When the line 
is blending, its contribution is "arbitrarily" situated and 
has only a little effect. In other words, the blending lines 
only contribute to the noise, perhaps modifying its statisti- 
cal distribution. From this point, noise is reduced when the 
number of spectral lines increases. 

In the above mentioned works, some 200 spectral lines 
were added in order to yield significant detections of 
Zeeman signatures. This straightforward method of line ad- 
dition was later improved by Semel (1995) in which a list of 
spectral lines of interest was created. These spectral lines 
were labelled by wavelength, Ai, equivalent width, u>i, and 
effective Lande factor, ^j. In the next step, each spectral line 
was represented by a Dirac function Wi ■ gi ■ 5{X — Xi) and 
transformed, as explained in Semel (1989), to the Doppler 
coordinate X, with dX = cdX/X, being c the velocity of 
light. Using u}'^ = WiQi, the convolution of the observed 
circular polarisation spectrum V{X) and the line list results 
in a MZS given by: 



Z-LA{Xj) = 



,LA 



(1) 



where Xi is the Doppler coordinate of the centre of each 
spectral line and Xj is the Doppler coordinate axis. The 
symbol Z\,\{xj) stands for the MZS due to line addition. 
As a rule, the convolution by means of a Fourier transform 
is found to be quite efficient. The line addition technique 
has also been applied with unit weights, i.e., lo]"^ = 1, with 
very similar results to the more refined version presented 
above: 



N 



Xi), 



(2) 



where SLA stands for simplified line addition. 

Although initially it was only used for the circular state 
of polarisation, since the MZS resulting from line addition 
is model-independent, it can be applied cither to Stokes /, 
Q, U or V . The reason is that no assumption have been 



made in order to compute the MZS, being just the brute 
force line addition. Therefore, the line addition is a very 
powerful and robust method to detect Zeeman signatures 
in stellar polarised spectra where the signal-to-noise ratio 
per individual spectral line is very poor. 

2.2. Least-squares Deconvolution 

The LSD technique was introduced for the detection of 
Zeeman signatures in stellar polarised spectra by Donati et 
al. (1997). They postulated that, in a given observed Stokes 

spectrum, the shape of the local Zeeman signature in circu- 
lar polarisation is identical for all spectral lines. This means 
that individual Stokes V profiles are assumed to correspond 
to a common basic Zeeman signature that is modulated by 
a proportionality factor. This coefficient is then fixed in the 
frame of the weak field approximation. In other words, they 
assume that the magnetic field is weak enough so that the 
Zeeman splitting is negligible as compared to the Doppler 
width of the spectral line. As a consequence, the radiative 
transfer equation can be solved analytically and it can be 
demonstrated that the Stokes V profile is proportional to 
the derivative of the intensity profile (Sears 1913; Landi 
degl'Innocenti & Landolfi 2004). Following this approxima- 
tion, Donati et al. (1997) foimd that the proportionality 
factor, ojf^^ = Xidifji, i.e., just the product of the cen- 
tral wavc;k;ngth of the transition. A;, the central depth of 
the line, di, and the effective Lande factor, gi. Donati et 
al. (1997) also normalised the weights to a mean value of 
500 nm. Based on all these assumptions, a particular MZS 
is extracted by means of a least squares method from the 
observed spectrum. 

We note in passing that u}'^^ and cof^, although appar- 
ently different, have a very similar behaviour. If the spec- 
tral line is assumed to be of Gaussian shape whose width 
is dominated by thermal effects, it is possible to show that 



,LA 



oc wf^^ in general. 



In more detail, the MZS obtained through the least- 
squares deconvolution method of Donati et al. (1997) is 
calculated solving the following weighted linear regression 
problem: given N lines with Stokes V profiles V{xj — Xi) 
sampled at the A^^ Doppler coordinates xj, find the pro- 
file Z{xj) that better approximates the observations in the 
least-squares sense given the fixed values of {w]^^'^}. In 
other words, it can be seen as the problem of fitting Nx 
independent straight lines whose slopes are given by the 
values of Z{xj) and the abscissae are the values of the 
weig hts {wP^}. The solut ion is obtained, imder the case 
of Gaussian noise, by minimising the following merit func- 
tions: 



X, 



: [v{x, - Xi) . 



Z{x,)\ 



(3) 



where the weights are usually chosen to be Wi = l/a{xj — 
Xi)^, the inverse variance of the expected noise at each 
velocity point. The solution to the previous problem is: 



Zi^snixj) = 



,LSD 



V{xj - X, 



EiWi{u;f^^)^ 



(4) 



In combination with maximum entropy codes (Brown et 

al. 1991, Donati & Brown 1997), LSD-based Stokes / and 
V profiles have been used for the production of magnetic 
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Fig. 1. The left panel shows the probability distribution of the normalised LSD weights, ojf^^ / J^ii'^f^^)'^ ■ The vertical 
lines indicate the position of for each temperature, where N is the number of lines included in the histogram. The 
right panel presents the probability distribution function of the normalised simplified line addition weights, 



surface maps of quite a large number of stars with great 
success. 

The main disadvantage of the LSD method relies on 
the assumptions made to obtain the MZS. It is not true 
in general that, locally, all spectral lines have a common 
shape (see for instance Kurucz 1993). Moreover, the weak 
field approximation used to obtain the proportionality co- 
efficients states that Stokes Q and U are zero to first order. 
This would make this technique not applicable to recover 
the MZS of the linear states of polarisation. 

2.3. Statistical properties of weights 

A fundamental point to clarify is why LSD works even 
though it is based on assumptions that may be difficult to 
admit. In order to shed some light on this issue, it turns out 
essential to analyse the statistical properties of the weights. 
The left panel of Fig. [T] presents the probability distribu- 
tion function of the normalised weights i^f^^ / ^iii^f^^)"^ 
for many lines in several line-lists obtained for different stel- 
lar temperatures. We have verified that the distribution is 
quite insensitive to the value of the surface gravity and the 
only strong dependency is seen with respect to the effec- 
tive temperature. The histograms have been calculated for 
lines above 300 nm and whose calculated line depth is larger 
than 0.4, which is the typical selection criterium used in the 
standard application of LSD. The vertical lines in Fig.[T]are 
the value of the inverse of the number of lines included in 
the histograms, in other words, the value of the weights 
in the simplest version of line addition given by Eq. 
Therefore, the LSD weights can be considered, to a good 
approximation, to be positive random variables whose dis- 
tribution is close to log-normal (which can be not badly rep- 
resented by a Gaussian distribution), centred around 1/N 
and with a certain dispersion. The log-normal shape arises 
because the product of several independent random quan- 
tities converge towards such a distribution, similar to the 
convergence towards a Gaussian distribution when adding 
many independent quantities. In our case, the weights are 
the product of three quantities that do not present a clear 
dependency (surely not the effective Lande factor and the 
rest of quantities). 



The right panel of Fig. [T] presents a similar statistical 
analysis but for line addition. In our study we use standard 
line-lists applied to LSD in which the equivalent width is 
not present. For this reason, we plot the weights for the sim- 
plified case in which the weights are just the effective Lande 
factors of all included lines, so that we plot the histogram 
of gi/J2i9i- The vertical lines again indicate the value of 
1/7V for each temperature. The weights can again be repre- 
sented approximately as Gaussian random variables centred 
around 1/A^ and with a certain dispersion. 

Assuming wavelength-independent Gaussian noise in 
the observations characterised by a variance cr^, the ob- 
served Stokes V{xj — Xi) profiles can be considered to be 
random variables characterised by a normal distribution 
with a mean value of ^i{xj) and variance cr^. In our case, 
fJ-{xj) is the hidden Zeeman signature at Doppler coordinate 
Xj that we want to recover. According to the results pre- 
sented in Fig.[Tl the MZS of Eqs. H]) and ([4]) can be approx- 
imated to be the addition of iV products of two Gaussian 
random variables. It is known that the probability distri- 
bution of the product of two Gaussian random variables X 
and Y is, in general, very complex to calculate. However, it 
is possible to obtain the mean and variances of such product 
through the calculation of the moment generating function 
(Craig 1936): 

E{XY) = fixf^Y 

V{XY) = fij,a^ + p.lcj\ + a\al, (5) 

where the p.{x,Y} represent the means of the two vari- 
ables while the a'^-^ y] their corresponding variances. 
In our case, X corresponds to the normalised weights, so 
that /ix = 1/N and (t\ = crj, with cr^ the standard devia- 
tion of the weights. Likewise, Y corresponds to the observed 
Stokes V profiles, so that /iy = ^^{xj) and cry = cr^. For cal- 
culating the MZSs, we have to add many such products. At 
this point, following the central limit theorem, the resulting 
probability distribution function can be calculated as the 
convolution of N Gaussian distributions whose mean and 
variance are given by Eq. ([5]). The resulting mean is just N 
times the mean of each Gaussian while the resulting vari- 



500 G 



Meir Semel et al. : Multiline Zeeman Signatures 
1500 G 



2500 G 



500 1OO0 1500 2000 2500 

wavelength (pixel) 



1OO0 1500 2000 2500 

wavelength (pixel) 



1OO0 1500 2000 2500 

wavelength [pixel] 






wavelength (pixel) 



wavelength (pixel! 



wavelength (pixel) 



Fig. 2. Top panels: spectral synthesis of the 9 selected lines of the multiplet 816 of Fei (see Tab. [TJ line 1 to 9 from left 
to right) having a spectral resolution of 3x10^ and a sampling of 100m s^^ (corresponding to steps of about 2.1mA in 
wavelength). From left to right the magnetic field strength is B = 500, 1500 and 2500 G. The inclination of the magnetic 
field vector with respect to the hne of sight is the same for all the columns, having a value of 50°. Bottom panels: The 
MZS resulting from the addition of the nine spectral lines. In all panels, the upper curves in solid line show Stokes / 
profiles while the lower curves in dashed line represent Stokes V. The thin lines are the MZS having a resolution of 
3x10^. The thick lines represent the MZS smoothed to reduce the spectral resolution to 75000. Additionally, the bottom 
panels also show the derivative of the smoothed low resolution intensity profile in dotted line. 



ance is N times the variance of each individual Gaussian. 
Therefore: 



with variance 



(6) 



(7) 



The signal-to-noise ratio increases roughly linearly with the 
number of added lines but reaches a maximum value at: 



^ rn ax — 



1 



{S/Nf 



(8) 



real 



where (S'/A^)rcai = tJ-{xj)/a is the real signal-to-noise ratio, 
i.e. the ratio between the amplitude of the hidden Zeeman 
signature and the noise. Equation ([8|) demonstrates that 
TVniax represents the maximum number of lines one should 
add in order to have a net increase in the signal-to-noise ra- 
tio. If more lines are added, the signal-to-noise ratio starts 
to degrade because of the intrinsic dispersion in the weights. 
The larger the ratio ^{xj)/a, the smaller the number of 
lines one should add until reaching the regime of S/N degra- 
dation. Likewise, the larger the dispersion of the weights, 
the smaller the number of lines one should add. 

The previous results suggest the surprising possibility 
of using purely random weights in the application of line 
addition or LSD. The only requisite is that the weights are 
positive definite with a fairly symmetric distribution peak- 
ing close to 1/A^ with a sufficiently small dispersion not to 
dominate in Eq. ([7]). Such a distribution should lead to re- 
sults similar to those presented in Eqs. ([6]) and ([7|). Also 



surprising is the fact that the analysis we have carried out 
suggests that the best solution is just to add the lines with- 
out any weight. In this case, no saturation due to the pres- 
ence of dispersion in the weights appears. This is a conse- 
quence of the fact that, since the observations are assumed 
to be described by the Gaussian distribution N{fi{xj), cr^), 
the maximum-likelihood estimation of its mean is exactly 
given by Eq. 

Summarising, one can consider that LSD and complex 
line addition schemes work because they give a response 
similar to that of the simplest line addition method. The 
fundamental reason is that spectral lines constitute a well- 
behaved statistical sample. The line-to-line differences are 
small enough (for instance, the dispersion in the value of 
the LSD or LA weights) so that the plain addition of lines 
makes the signal appear as soon as one adds a sufficient 
number of lines. Although the assumptions on which LSD 
or any method are not really fulfilled, the addition of lines 
dilutes any possible difference in their behaviour. 



3. The magnetic sensitivity of the MZS 

The MZS are very efficient detectors of the magnetic ac- 
tivity of the star by means of the Zeeman effect. However, 
since our main interest is to understand the magnetism of 
the stars, we need to infer the magnetic field vector of the 
star from the information contained in the polarised spectra 
or from the MZS. In cools stars the signal to noise ratio per 
individual spectral line is usually so poor that it is impos- 
sible to make reliable measurements of the magnetic field. 
Consequently, we are obliged to use the MZS as indicators 
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Table 1. The list of spectral lines of multiple! 816 of neu- 
tral iron. The relative strengths correspond to LS coupling 
as given by Allen (2000). The ratio rj of line to continuum 
absorption is obtained from the relative strengths multi- 
plied by a factor of 2; the latter is chosen arbitrarily. 



of the magnetic activity in the stars. At this point a ques- 
tion arises: does the MZS contain the information of the 
magnetism of the star? It is evident that, in the case of 
the line addition, the MZS is not a physical spectral line. 
Then, does it behave like a spectral line and contain any 
information about the magnetic field?. This section will be 
devoted to study to what extent it is possible to use the 
MZS to reliably infer the magnetic field in stars. 

We simulate the observed spectral lines from a star as- 
suming a Milne-Eddington atmosphere. Moreover, we sup- 
pose that there is only one magnetic field vector in the stel- 
lar surface that is responsible for the Zeeman signatures 
in the spectra. Although more sophisticated scenarios and 
solutions to the polarised radiative transfer equation exist, 
we are not interested in a realistic simulation of the stellar 
spectra but in showing that, once the information of the 
spectra is encoded in the MZS, we are capable of retrieving 
it. Using the M-E model facilitates the analysis and does 
not require heavy calculations. We use the code Diagonal as 
described in Lopez Ariste and Semel (1999) with only one 
layer. This code can arbitrarily cope with a large number 
of layers if needed. For the sample M-E model atmosphere 
we choose the value of the gradient of the source function 
to be /? = 10 and the ratio between the line and continuum 
absorption, r?, to be equal to twice the relative strength of 
the components of the multiplet as given by Allen (2000). 
Doppler broadening is set to 21 mA as in sunspots, and the 
reduced damping is taken to be 0.01 Q 

We synthesise the Stokes vector of only a few spectral 
lines with a spectral resolution close to 3x10^. We select 
the multiplet 816 of Fei (listed in Tab. [I]) to ensure the va- 
lidity of LS coupling and the easy determination of Zeeman 
patterns, relative line strengths, etc. We then add up the 
Stokes parameters of all the spectral lines using unit weights 
[therefore, using the simplified line addition of Eq. and 
subsequently reduce the spectral resolution to stellar con- 
ditions, say 75000, obtaining the MZS. 



^ Saturation has always been an important issue in solar and 
stellar magnetic field measurements. Here it is included in the 
procedure and we have to use the same procedure when we get 
to the inversion. 



3.1. MZS in circular polarisation 

Figure [2] shows the synthetic Stokes / and V together with 
the MZS for the selected Fei lines for three different values 
of the magnetic field strength. We keep a high spectral res- 
olution of 100 m/s, i.e. a resolution of 3 millions, and then 
study the effect on the MZS of smoothing, bringing the res- 
olution down to 75000, which corresponds to sampling with 
a step of about 4 km s~^. 

With very high spectral resolution both Stokes / and 
Stokes V profiles change morphology and details clearly 
depending on the field strength. In the case of 500 G we do 
not detect a clear magnetic splitting in the Stokes / profiles. 
Also, the Stokes V profiles are all very similar, the only dif- 
ference being a scale factor. This is a consequence of the 
fact that for 500 G, all spectral lines are in the weak field 
regime of the Zeeman effect, where the magnetic splitting is 
smaller than the Doppler width of the line. However, when 
the magnetic field increases, the Zeeman splitting becomes 
evident in the intensity profiles and the circular polarisation 
profiles abandon the shape-similarity regime. In fact, the 
anomalous dispersion terms are now important and show 
the typical reversal of polarities in the centre of the Stokes 
V profile. This has no counterpart in the weak-field approx- 
imation since anomalous dispersion is a high-order effect. 
Although the spectral lines are now out of the weak field 
regime, they have not still reached the strong field regime 
(where the Zeeman splitting is very large as compared to 
the Doppler width of the line) . 

Considering the MZS, in the weak field case (500 G), 
we see that, even in the high spectral resolution case, the 
shapes of the intensity and circular polarisation MZS also 
resemble the individual Stokes / and Stokes V profiles, re- 
spectively. When lowering the spectral resolution to the 
typical stellar one, the profiles are also similar in shape 
but the amplitude has decreased and the profiles are wider. 
Note that the MZS in circular polarisation matches almost 
perfectly with the derivative of the intensity profile. The 
reason is simple to understand. Since all the lines are in 
the weak-field regime and the derivative is a linear opera- 
tion, the addition of derivatives equals the derivative of the 
addition of lines. This confirms that in this case the weak 
field regime is a good approximation. 

When we leave the weak field regime, the MZS at high 
spectral resolution shows some distortions that are due to 
the different shapes of all the individual Stokes V profiles 
(every line departs from the weak-field regime at different 
values of the magnetic field strength). When we decrease 
the spectral resolution to the stellar case, in the case of 
1500 G we see that the MZS for circular polarisation is 
close to the derivative of the intensity MZS profile. This 
means that, even if the individual spectral lines are not 
in the weak field regime, the MZS for circular polarisation 
can be in this regime when lowering the spectral resolu- 
tion. The reason behind this behaviour is that the width of 
the lines are increased when lowering the spectral resolu- 
tion, thus extending the validity of the weak-field regime. 
Consequently, the weak field regime applies for larger mag- 
netic fields in the circular polarisation MZS than in the 
individual spectral lines (typically up to 600-700 G in the 
visible part of the spectrum). However, this does not occur 
in the 2500 G case. The low resolution circular polarisation 
MZS already differs substantially from the derivative of the 
intensity MZS. 
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Fig. 3. The MZS for a magnetic field strength of 1000 G (left panels) and 3000 G (right panels). Inside each panel, the 
intensity profile are shown at the top, while the polarisation profiles are shown at the bottom, with Stokes V shown in 
dashed lines and Stokes Q in solid lines. High spectral resolution profiles are shown in thin lines, while the smeared low 
spectral resolution profiles are indicated with thick lines. 



The MZS constructed with line addition are not sim- 
ple spectral lines in the sense that they have no associated 
wavelength and are not characterised by specific atomic pa- 
rameters such as the excitation potential. Moreover, in the 
more general case they do not even correspond to a partic- 
ular chemical element, and no specific spectroscopic term is 
attached to them. However, the Zeeman signatures do not 
disappear in the MZS. We have seen that they are sensi- 
tive to the magnetic field, having also a weak field regime 
where the amplitude of the circular polarisation MZS scales 
with the derivative of the intensity MZS. For this reason, 
they can be used to extract information about the magnetic 
field, as we show in Sg] 



3.2. MZS in linear polarisation 

As a rule of thumb, the Zeeman linear polarisation is at 
least one order of magnitude lower than the circular one 
(e.g. Carroll et al. 2007). Now, if our line of sight (hence- 
forth LOS) is not a preferred direction for the stellar mag- 
netic field vector and therefore its three components (the 
longitudinal one and the two perpendicular to the LOS) 
have all equal probabilities, on average the transverse com- 
ponent is times stronger than the LOS component. 
From solar physics we know that kG fields are not only 
found in sunspots. However, can we assume that kG fields 
are common also in solar type stars? If yes, why does the 
linear polarisation escape observation? In the following, we 
try to answer this question and we also suggest a remedy. 



Let us denote ip the angle between the magnetic field 
vector and the LOS and take it between 0° and 90°, say pos- 
itive LOS component. If all directions have the same prob- 
ability, the average inclination angle is tp = Irad ~ 57.3°. 
Likewise, the median angle is f/'Median = 60°. In other 
words, there are as many orientations between 0° and 60° 
as between 60° and 90°. For the case of negative J^OS lon- 
gitudinal fields, 90° < -0 < 180°, the results for tp and for 
ip Median Stve dcduccd similarly. In conclusion, ip is likely to 
be nearly 60°. 

Figure [3] shows the MZS for the Stokes /, Q and V pa- 
rameters computed with a magnetic field strength of 1 kG 
and 3kG and an inclination ip = 60°. Note that, for high 
spectral resolutions, the MZS for Stokes Q is approximately 
0.7 times the MZS for Stokes V. Stokes Q changes with 
wavelength twice as fast as Stokes V, so that for current 
high spectral resolutions, say 75000, Stokes Q shrinks much 
faster than Stokes V in the convolution process. Moreover, 
Stokes Q changes sign when the azimuth changes by 90°, 
while Stokes V changes sign when tjj changes by 180°. Both 
effects can be reduced by increasing the spectral resolution, 
to values of the order of 120000. This means that observing 
at high spectral resolutions. Stokes Q is preserved consid- 
erably and, more interestingly, the spatial resolution of the 
stellar surface is improved as well. 

4. Recovering the magnetic field from the MZS 

Our purpose is to assess to what extent the new MZS still 
contains all the information on the magnetic field that was 
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Fig. 4. Top; A test of the centre of gravity method. The 
CoG shifts are calculated over a database of 30 values of 
magnetic fields in steps of 100 Gauss from to 3kG. The 
shifts are divided by the cosine of the inclination angle and 
by the adopted effective Lande factor to compare with the 
amplitude of the magnetic field (abscissa). The curves are 
shown for the following values of V', the inclination angles in 
degrees: (solid line), 20 (points), 40 (dashed), 60 (dashed 
- points) and 80 (stars). The value of g for the pseudo-line 
was chosen to be 1.59, For = the CoG shifts correspond 
exactly to the line of sight component of the magnetic field. 

present in the individual spectral lines in our sample. To 
this end, we examine the invariants that we have already 
found in solar magnetism, like the displacement of the cen- 
tre of gravity (CoG) that yields a good approximation to 
the longitudinal component of the magnetic field. It is im- 
portant to have in mind that the centre of gravity is a linear 
operator and commutes with the algebra of line addition 
and with smoothing. 

The centre of gravity method applies only to circular 
polarisation; it is well studied in solar magnetism and de- 
scribed in a number of papers of which at this point we 
mention only Rees et al. (1979) (and references therein). In 
the limiting case of optically thin layers in local thcrmody- 
namical equilibrium, the centre of gravity shift of the line 
profiles observed in circular polarisation is proportional to 
the longitudinal component of the magnetic field. This still 
holds true, even for line formation in optically thick layers, 
in the limiting case of weak magnetic fields. Usually, this 
simple method is still a fair approximation in the more gen- 
eral case, but there are a few exceptions - see Semel (1967) 
for theoretical demonstrations, experimental tests and some 
exceptions. For further discussion of the CoG method ap- 
plied to stars see also Stiff (1986) and Leone & Catanzaro 
(2004). 

We apply the centre of gravity method to the MZS of 
Stokes V (the ones showed in Fig. [2]). We compute the MZS 
for Stokes Y assuming a magnetic field strength going from 
to 3 kC having different inclinations t/' = 0, 20, 40, 60 and 
80°. Then, we infer the magnetic field following the cen- 
tre of gravity method. Figure [4] shows the plot of the de- 
duced magnetic field versus the theoretical one. The figure 
shows a clear linear trend showing that the MZS contains 
the information of the magnetic field that was present in 
the individual spectral lines. Since the centre of gravity is 



applicable both to the individual spectral lines and to the 
MZS of Stokes V ^ we conclude that the MZS encodes the 
information of the magnetism of the star in a similar way 
that the individual spectral lines do. In other words, the 
MZS constructed with the line addition technique behaves 
as a spectral line although not any specific atomic element 
is associated to it. However, it can be characterised by an 
effective lande factor and the extraction of the magnetic in- 
formation is possible. Figure |4] also shows an error close to 
10% of the inferred value which is intrinsic to the centre of 
gravity method. Our aim here is to show that the magnetic 
information is present in the MZS and not to do precise 
measurements of the magnetic field strength. In order to 
do these exact measurements we must use more sophisti- 
cated techniques. To have reliable values of the magnetic 
field, it is fundamental that the inversion procedure uses 
exactly the same operator that generates the MZS for the 
comparison with the observational data. This topic will be 
discussed in a forthcoming paper. 

5. Conclusions and discussion 

We have used the line addition technique to show how to 
extract a Zeeman signature for any of the Stokes parame- 
ters. These technique, apart from being applicable to any 
state of polarisation, is also model independent. In other 
words, we do not need to do any approximation nor any 
assumption on the magnetic field or the stellar atmosphere 
to extract the Zeeman signature from the spectra. The line 
addition, as simple as the brute force addition of many spec- 
tral lines is a powerful tool to detect the magnetic field in 
cool stars, where the Zeeman signature per spectral line is 
very poor. 

Moreover, the MZS created by the line addition tech- 
nique also contains the magnetic information of the star. 
Therefore, it is possible to infer the magnetic (and the ther- 
modynamical) properties of the stellar atmosphere using 
the appropriates inversion procedures. We have shown that 
the magnetic field strength can be inferred from the MZS 
of Stokes V by using a procedure as simple as the cen- 
tre of gravity. However, we must point out that, in order 
to have more reliable measurements of the magnetic field 
vector, we need to apply more sophisticated inversion pro- 
cedures. Such procedures should ideally take into account 
in the process the same operator used to create the MZS 
when comparing with the observations. 

Linear polarisation is often neglected in stellar polarime- 
try under the assumption that it will present amplitudes too 
small to be observed. In this paper, we have demonstrated 
that, at appropriate spectral resolutions, the linear polari- 
sation due to the Zeeman effect produces sensible signatures 
useful for diagnostics of the stellar magnetic fields. With a 
resolution above 10^, we can probably detect significant sig- 
nals in linear polarisation in the spectra of solar type stars. 
Such a resolution has been achieved in spectroscopic tests 
with the SemPol instrument working together with UCLES 
at the A AT with the 79G/mm gratin£|. The detection of 
linear polarisation in the spectrum of a solar type star has 
been performed in 2004 (see Semel et al, 2006). The linear 
signal observed was indeed four time less than the circular 

^ See the discussion by by Lopez Ariste & Semel 
(1999) on the spectral resolution of UCLES in 
http: //www . aao . gov . au/local/www/UCLES/ cookbook/report2 . 
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one, but still significant. There is definitely an interest to 
proceed in this direction. 

We have also demonstrated why line addition and least- 
squares dcconvolution give very similar MZS. The funda- 
mental reason is that the weights used in both cases have 
well-defined statistical properties: their probability distri- 
bution is positive definite, centred at the inverse of the 
number of lines and not too broad. These properties induce 
that weighted line addition and LSD converge to the sim- 
plest possible line addition which, since noise is assumed 
to be of Gaussian type, turns out to be the maximum- 
likelihood estimation of the MZS. This opens the curious 
possibility of applying LSD with random weights provided 
they fulfil the mentioned properties. 

In this paper we have presented a very simple technique 
to build the MZS. However, more powerful techniques are 
being developed in the frame of PCA. The main idea is to 
build a more sophisticated operator than the simple Dirac 
functions used for the line addition. The operator using 
PCA will be related to the eigenvectors of a data base 
containing theoretical stellar spectra synthesised in sev- 
eral magnetic and thermodynamic configurations. Ramirez 
Velez et al. (2009, submitted) explains this technique in 
detail and apply recent solar inversion procedures to the 
synthesised MZS using PCA. Other applications of PCA to 
obtain the MZS has been presented by Martinez Gonzalez 
et al. (2008) and by Carroll et al. (2007). Apphcation of 
PCA has been proven very effective in solar magnetome- 
try and, analogously to the solar database will be 

created and the PCA eigenvectors derived, applying the 
PCA-ZDI detection to just one magnetic point on the star. 
In later stages of the work, we shall discuss how to apply 
our approach to some few points on the stellar disc well sep- 
arated by significant Doppler effects and, finally, we shall 
treat the more realistic case of continuous distributions of 
fields over the stellar surface. Here, the orthogonality of 
the eigenvectors is not conserved for adjacent stellar points 
(that is for small differences in the Doppler shifts). 
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